Supersolitons: Solitonic excitations in atomic soliton chains 
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We show that, by appropriately tuning physically relevant interactions in two- component nonlin- 
ear Schrodinger equations, it is possible to achieve a regime with particle-like solitonic collisions. 
This allows us to construct an analogue of the Newton's cradle and also to create localized collective 
excitations in solitary-wave chains which are quasi- integr able solitons, i.e. supersolitons. We give a 
physical explanation of the phenomenon, support it with a perturbative analysis, and confirm our 
predictions by direct simulations. 
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Introduction and model.- One of the most successful 
concepts of nonlinear science with applications to a great 
variety of physical contexts is that of solitons, i.e., self- 
localized nonlinear waves, sustained by the balance be- 
tween dispersion and nonlinearity. Many types of solitons 
have been studied, from the classical examples found in 
integrable models, such as the Korteweg - de Vries, sine- 
Gordon, Toda-lattice (TL), nonlinear Schrodinger, and 
other celebrated equations, and extending into the realm 
of realistic non-integrable nonlinear- wave models. 

Solitons are usually expected to be robust against col- 
lisions, which is a trademark feature of integrable equa- 
tions. A lot of activity has been directed at the study of 
soliton collisions and interactions in non-integrable sys- 
tems. Recent advances include the analysis of chaotic 
scattering [1], the formation of soliton bound states and 
soliton clusters [2] ^nd the studies of soliton collisions in 
vector systems [ll, Ij, ^ , to name just a few. 

While it is customary to speak of solitons as elastically 
colliding quasi-particles, most solitons, specifically in in- 
tegrable systems, pass through each other, thus clearly 
featuring their wave nature. On the other hand, elastic 
collisions between classical particles lead to momentum 
exchange between them, and rebound, due to the non- 
penetrability of classical particles. 

In this paper we discuss a particular soliton colli- 
sion scenario of physical relevance, where truly elastic 
particle-like soliton collisions can be achieved. We will 
show how this can be used to build a vector-soliton ver- 
sion of the Newton's cradle, and to build supersolitons, 
i.e. collective soliton-like excitations in arrays of solitary 
waves, leading to a remarkable conjunction of emergent 
phenomena: the former one representing the formation 
of robust soliton trains, and the latter effect implying 
the emergence of an effectively quasi-discrete soliton at 
a higher level of organization. 

The basic model which allows to implement the above- 
mentioned effects is based on the two-component (vecto- 



rial) nonlinear Schrodinger equation (NLSE), that arises 
in sundry contexts [5]. Its normalized form is 
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An important physical realization of this model is a mul- 
ticomponent Bose-Einstein condensate (BEG), where Uj 
are wave functions of two atomic states under the action 
of a strong transverse trap with frequency z/_l p. The 
variables x and t are measured, respectively, in units of 
ao = \/Iilmv^ and l/z^±, and gij = 2aij/ao, with aij the 
respective s— wave scattering lengths. The normalization 
integral for Uj gives the number of atoms in the respective 



species, /_ 



^(fx = JVj. Solitons in BEGs have been 



created experimentally [7| and their interactions studied 
theoretically in many papers (see e.g. jslJoUlOt). 

Particle-like elastic collisions and the solitonic New- 
ton's cradle.- We will consider Eqs. ([1]) with intra- 
component attraction (^11,^22 < 0) and inter-component 
repulsion (^12,^21 > 0). To fix ideas we will choose 
gii = ^22 = —912 = — ^2i(= 1 without loss of generality) 
for which case the coupled NLSE are not integrable |ll| . 
In this situation the solitons, that may be formed in both 
components independently, interact incoherently with a 
repulsive force. The basic physical feature underlying our 
analysis is that the dynamics of those solitons may be 
similar to that of elastic beads. We will explore the cases 
of harmonic longitudinal confinement V{x) = Q'^x'^ /2, 
and ring-shaped configurations [12]. 

Let us consider soliton trains built as follows: 

ui{x,0) = 2_. (— l)^sech(x — ^n)exp (zx'Un) (2a) 
^2(^,0) = ^ (-l)'^sech(x - Cn)exp(zxK;n)(?b) 
with alternation of the soliton species in the train, i.e., 

•••^n-l < Cn-1 < £,n < Cn < ^n+1 < Cn+1 < •••, (3) 



Vn and Wn being initial velocities of the solitons. In Fig. 
[ll where the trap is absent, panel (a) displays a single 
collision event {N =1). It is noteworthy that, because 
of the repulsive inter-component interaction, the incident 
soliton (in field ui) transfers all of its momentum to the 
initially static soliton (in 2x2), in full compliance with the 
behavior of elastic particles and contrary to the typical 
behavior of nontopological solitons in integrable systems. 
The dynamics of a train of eight alternating solitons in a 
ring configuration demonstrates the periodic transfer of 
the momentum through the train, see Fig[T](b). 




FIG. 1: [Color online] Pseudocolor plots of \ui{x,t)f (blue) 
and \u2(x,t)\'^ (red) show the time evolution of initially well 
separated soliton trains built as per Eq. (|2]). The horizontal 
(vertical) axis of each plot is spatial adimensional variable x 
(adimensional time t). (a) Collision of two solitons (N = 1) 
for ^1 = -20, Ci = 0,vi = 0.4,^1 = 0, X G [-30,30] and 
t G [0, 100]. (b) Multisoliton collisions in a ring for A^ = 4, 
(j = -35 + 20j,^^- = -25 + 20J, x G [-40,40] and t G [0,250]. 
Initial velocities are all zero except for W3 = 0.5. 




FIG. 2: (Color online) Same as Fig. [T]for a Newton's cradle 
built up of five solitons with an external potential with Q^ = 
6 • 10"^ In all subplots tmax = 1000, x G [-40,40], and the 
initial separation between solitons is ^n — Cn =6. Different 
oscillation modes are excited by applying an initial velocity, of 
size 0.2, to: (a) the top-left soliton, (b) two ultimate solitons 
on both sides, with opposite velocities, and (c) all solitons. 

By adding the parabolic trapping potential we urge 
the solitons to oscillate around the equilibrium position. 
In the quantum interpretation of our model this setup 
would provide a way to construct a quantum Newton^s 
cradle with atomic solitons, as illustrated in Fig. [2l Un- 
like other settings explored in BEG [l3], the cradle con- 
figuration does not require a lattice potential to create 



effective particles, which are here created by purely non- 
linear interactions. 

The Toda-lattice limit: super solitons - We now con- 
sider two alternating chains of solitons set along a ring. 
Within the effective-particle approach, we assume that all 
solitons, which have identical amplitudes in each compo- 
nent, T] (for ui) and (for 1^2), behave like rigid particles 
and thus do not suffer conspicuous deformation, i.e., each 
soliton may be approximated by 

u^^\x,t) = 7?sech[7/(x-^n)]e'^^''+^^('^'"^')^',(4a) 



U2 {x^t) = Osech.[0 {x — (n)]e 



iC^x+^J{e^-C?,)dt 



(4b) 



with the initial positions arranged as per Eq. (|3]). 

A straightforward analysis based on the perturbation 
theory for solitons [14] yields the following system of 
equations of motion for the soliton coordinates: 



Cn 



Cn 



-n^Cn + 8r] [r?2 (e- 
+0'[ 
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These equations are derived under the assump- 
tion that adjacent solitons are well separated, i.e., 
{V^ ^) {in - Cn-i) , (^, 0) {Cn - ^n) > 1, although a strong 
inequality is not really necessary here. Similar ideas 
have been used to derive equations for the interaction of 
other elementary nonlinear structures, which gives rise 
to different equations at a higher level of organization 
0, HO, [l5|]. Equations ^ with 1] = reduce to the 
so-called diatomic TL, which is not integrable, although 
some solutions are known. 

With T] = and 1] = in Eqs. ([5]) and defining q2n(f) = 
'^V(.n{t),q2n-\-i{t) = 2r]Cr,it) and a = 32/^^, we arrive at 
the integrable TL model [l6|, 



-iqn-qn-l) 



-{qn+l-Qn) 



(6) 



This model describes a dynamical lattice with the ex- 
ponential potential of the nearest-neighbor interaction. 
However, potentials of the interaction between adjacent 
atoms in real condensed-matter systems are never expo- 
nential, being closer to those of nonlinear anharmonic 
oscillators. This is why the only experimental realization 
of the integrable TL was realized in electric transmission 
lines Lis,], that may be readily designed in exact corre- 
spondence to Eqs. dn]). Our analysis suggests a possibil- 
ity to create Toda solitons, of both mono- and diatomic 
types, as excitations in interwoven arrays of multicom- 
ponent NLSE solitary waves. 

We name these excitations supersolitons since they oc- 
cur on top of an array of "elementary" solitons, and are 



expected to be as robust as solitons in integrable mod- 
els. The same name was previously applied to solitons in 
supersymmetric models [17], and, in a completely dif- 
ferent context, to localized topological collective exci- 
tations in chains of fluxons trapped in periodically in- 
homogeneous Josephson junctions and in layered super- 
conducting structures [l5|. The appearance of TL su- 
persolitons represents a remarkable phenomenon at a 
higher-organization level, using, as building blocks, soli- 
tary waves of the multicomponent NLSE, i.e. a strongly 
nonintegrable model. 

In the monoatomic lattice, Eq. (|6|) has an obvious equi- 
librium solution with Qn = L/ (2A^), where TV is the num- 
ber of solitons in each subchain, and L the total length 
of the system. For small perturbations with frequency 
uj and wavenumber k around this configuration, the dis- 
persion relation is cj^ = 1287/'^e-^^/^ sin^(fc/2). With 
respect to the quantization imposed by the boundary 
conditions for the ring-shaped soliton chain, k = Trm/N^ 
m = 0, ±1, ±2, ... , this yields 

\ujm\ = ^V2rf exp {-r]L/2N) |sin(7rm/(27V))| . (7) 

If a wave in the lattice is excited by kicking one soliton 
and lending it velocity v, the wave will hit solitons with 
period T = L/(2Nv)^ which corresponds to an effective 
excitation frequency c^exc = 27r/T = AiiNv/ L. Thus, res- 
onant excitations may be expected under the condition 
PuOe^c = Ql^ml, or, in other words, at values of the kick 
velocity belonging to the following resonant spectrum^ 



XP,Q) 



Q 2V2Lrf 



ttN 



-'^^^^^\sm{7rm/{2N))\, (8) 



where integers Q and P stand for the order of the reso- 
nance and subresonance {P = Q = 1 correspond to the 
fundamental resonance). Another interpretation of this 
resonance condition (cf. Ref. [19]) is that the kick veloc- 
ity coincides with the phase velocity of linear waves. 

Numerical studies of supersolitons - To verify our pre- 
dictions based on Eq. (|6]), we have performed numerical 
simulations of Eq. ([T]). First, in Fig. [3] we have gen- 
erated a single supersoliton by kicking one of the most 
external solitons in one of the components. Since this 
excitation does not correspond exactly to a supersoliton 
we also obtain a small ammount of radiation which is 
seen as small remnant oscillations of the individual soli- 
tary waves. Apart from this efect due to the excitation 
procedure, the propagation of the supersoliton is perfect 
as seen both in the amplitude [Fig [31(a)] and pseudocolor 
[Fig. [3jb)] plots. Another effect seen in Fig. [31(a) and 
not considered in our model is the small compression of 
the individual solitary waves when they are hit by the 
supersoliton (the model assumes equal amplitude indi- 
vidual solitary waves). However this small effect does 
not affect our conclusions and can be minimized by con- 
sidering smaller energy collisions (i.e. incident speeds). 



Fig. [H shows the collisional behavior for head on colli- 
sions of equal speed solitons [Fig. [l^a)] and the overtak- 
ing of a slow supersoliton by a faster one [Fig. [41(b)] . In 
both cases the supersolitonic excitations behave as true 
solitons, what it is justified by the integr ability of our 
simple model given by Eqs. (|6|). We want to emphasize 
again that these behaviors, typical of integrable systems, 
arise on top of a strongly nonintegrable model. 




FIG. 3: (Color online) Excitation of a supersoliton on a soli- 
ton chain given by Eq. ([2][3]) with N = 24 localized solitons 
in each component, by kicking the j — 22soliton with veloc- 
ity W22 = —0.5. Individual solitons have initial amplitudes 
?7 = ^ = 1 and the intersoliton distance is five units lead- 
ing to a total ring length of 240. Left: Amplitude plot of 
|^(x,t)|^ = |i^i(x,t)p + \u2{x,t)\'^. Right: Pseudocolor plot 
showing |i^i(x,t)|^ (yellow) and \u2{x,t)\'^ (red). The range 
of adimensional times spanned in both plots is t G [0, 250]. 
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FIG. 4: (Color online) Examples of supersoliton collisions. 
The initial configuration is as in Fig. [3]but now two supersoli- 
tons are excited. Left: Head-on collision induced with excita- 
tion parameters V3 — 0.5 and W22 — —0.5. The range of adi- 
mensional times spanned is t G [0,250]. Right: Soliton over- 
taking excited with initial velocities V23 = —0.6, wiq = —0.3. 
The range of adimensional times spanned is t G [0, 350]. 

Can scalar models support supersolitons?- Soliton col- 
lisions in the framework of scalar NLSEs have been stud- 
ied in various contexts , and equations similar to Eqs . 
([5]) have been derived using different approaches [SJ, [lO[ , 
leading to the so-called complex TL. Despite the formal 
similarities, the ensuing dynamics is not robust, and soli- 
tonic solutions turn out to be unstable because of the 



phase dependence of the interactions. An example is dis- 
played in Fig. EJa), which shows that the phase shifts 
induced by the initial kick velocity lead to an unstable 
dynamics of the single-component chain, whereas in its 
alternating two-component counterpart the system does 
not display any instability, as shown in Fig. [5jb); in par- 
ticular, the configuration shown in Fig. Efb) periodically 
recovers its shape. Thus, the vectorial system with in- 
coherent interactions is free of the instability of soliton 
chains in single-component models with coherent (phase- 
dependent) interactions [2^]. 







FIG. 5: (Color online) Time evolution for < t < 500 of 
matter-wave soliton trains with a small collective excitation 
in the case of one (a) and two (b) species, which are de- 
scribed, respectively, by the scalar and coupled NLSEs (black 
and white colors). Adjacent solitons are given opposed input 
velocities, x — ±0.1 in (a), and Vj = 0.1, Wj = —0.1 in (b). 
Other parameters are as in Fig. [S] 

Experimental realization.- The creation of TL super- 
solitons in BECs would depend on the use Feshbach res- 
onance techniques to get an atomic mixture with attrac- 
tive intra-species and repulsive inter-species interactions. 
Atomic mixtures with controllable interspecies interac- 
tions have already been reported in Ref. [21]. Our ini- 
tial state of alternating solitons may be created by the 
modulational instability and se greg ation from an initially 
stable two-component mixture [22| . 

Conclusions.- We have explored a physical model 
based on the vectorial NLSE, in which hard-particle-like 
(bouncing) elastic collisions between solitons belonging 
to different species are possible. These interactions allow 
building an analogue of the Newton's cradle using soli- 
tary waves, and to supersolitons in a chain of alternating 
solitons. The existence of these robust localized collec- 
tive excitations on top of arrays of nonintegrable solitons 
represents a remarkable emergent phenomenon. 
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